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ABSTRACT 

K2 space observations recently found that three super-Earths transit the nearby M dwarf K2-3. The apparent bright¬ 
ness and the small physical radius of their host star rank these planets amongst the most favourable for follow-up 
characterisations. The outer planet orbits close to the inner edge of the habitable zone and might become one of the 
first exoplanets searched for biomarkers using transmission spectroscopy. We used the HARPS velocimeter to measure 
the mass of the planets. The mass of planet b is 8.4 ±2.1 M 0 , while our determination of those planets c and d are 
affected by the stellar activity. With a density of 4.32lQ'7g g cm“®, planet b is probably mostly rocky, but it could 
contain up to 50% water. 
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1. Introduction 


Crossfield et al.| (2015 1 have recently reported from K2 
IHowell et al. 2014p observations that the star K2-3 


2MASS 11292037-0127173, EPIC 201367065) hosts a plan¬ 


etary system with three transiting super-Earths. The star 
is an inactive MO dwarf (Teff = 3900 ± 190 K, [Fe/H] = 
—0.32 ± 0.13 dex), and it is bright enough to be amenable 
to transit spectroscopy. With a radius of 1.5 R 0 , Planet d 
orbits close to the inner edge of the system’s habitable zone. 
The mass of the planets has not been measured to date. 


To constrain those masses, we monitored the radial ve¬ 
locity of K2-3 with the HARPS velocimeter (Mayor et al.| 


20031. We jointly analysed these new velocities and the 
K2 photometry through n-body integrations to obtain the 
mass of planets. Combined with the planet sizes, this con¬ 
strains their mean densities and bulk compositions and 


* Based on observations made with the HARPS instrument on 
the ESO 3.6 m telescope under the program ID 191-C0873 at 
Cerro La Silla (Chile). 


helps explore the rocky-gaseous transition in the super- 
Earth regime. 


2. Observations 

2.1. K2 light curve 

The K2 mission observed K2-3 during its Campaign 1 
(Summer 2014) in long cadence mod^ We downloaded the 
pixel data from the Mikulski Archive for Space Telescopes 
(MAST)[^ and used a modified ve rsion of the CoRoT im- 
agette pipeline ( Barros et al.|2014 ) to extract a light curve. 
The procedure first determines the circular synthetic aper¬ 
ture that maximises the photometric signal-to-noise ratio 
on the mean image. For each ima ge, it then computes a 
modified moment method centroid ( Stone|19^ I and recen¬ 
tres a heavily oversampled version of the image on the cen¬ 
troid before extracting the flux inside the pre-determined 
aperture. The degraded pointing stability of the K2 mis- 


^ http://archive.stsci.edu/k2/ 

^ http://archive.stsci.edu/kepler/data_search/search.php 
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sion couples with pixel sensitivity variations to introduce 
position-dependent systematics in the raw light curves. To 
correct for this flux d ependence with position , we u sed a 
procedure inspired by Vanderburg & Johnson (20141. The 
light curve was divided into seven equal duration segments, 
and for each segment we perfor med a ID decorrelation as 
described in Vanderburg ( |2014 1. We found that a 21-pixel 
photometric aperture results in the best photometric pre¬ 
cision of the final light curve with a 204 ppm mean RMS. 
The 80-day light curve contains eight transits of Planet b, 
four of Planet c, and two of Planet d. We only modelled the 
light curve around those transits, after normalising it with 
a parabolic fit to its out-of-transit part. To account for the 
29.4 minutes integration time of the long cadence data, we 
oversampled the model light curve by a factor of 20 and 
then binned it to the cadence of the data points. 


2.2. HARPS radial velocities 

We obtained differential radial velocities (RVs) with 
HARPS, the ESO velocimeter installed at the focus of the 
3.6m telescope at La Silla Observatory (Mayor et al. 2003; 
Pepe et al. 2003). We chose not to use simultaneous wave¬ 
length calibration and to instead rely on the <1 ms“^ sta¬ 
bility of the HARPS spectrograph, since we expected sig¬ 
nificantly larger photon noise errors. We tried to secure two 
1800 s observations per night and collected 66 spectra over 
a timespan of 103 days. 

For an optimal extraction of the velocity signal, we used 
all recorded spectra (already extracted and wavelength- 


G1 846^an M0.5V star with a 


calibrated by the online pipeline, Lovis & Pepe 20071 to 
build stellar templates. We shifted each spectrum by its 
barycentric correction and co-added those shifted spectra 
into higher S/N templates after carefully masking the tel¬ 
luric absorption lines. 

When using such a small set of relatively noisy spectra, 
we would bias the velocity if we included the spectrum anal¬ 
ysed for Doppler shift in the template that it is matched 
against, since their common noise pattern will contribute 
to the match. To avoid this bias, we built one template for 
each spectrum by co-adding all spectra but the spectrum 
itself. 

We then computed the RV shifts that minimise the 
difference between the individual spectra and their tem- 


plates (e 

.g. Zucker & Mazeh|2006 

Anglada-Escude & But- 

ler 

2012 

Astudillo-Detru et al.||2014|). The 66 resulting ve- 


median uncertainty of 2.9 ms“^ (compared with 7.1 ms“^ 
and 4.3 ms“^ for the velocities measured by the HARPS 
pipeline). Although the orbital periods are known from the 
photometry and do not need to be identified from the RVs 
alone. Fig. |A.l| shows their periodogram for illustration. 

2.3. Stellar activity 

The HARPS spectral range includes both Ha and the 
Can H&K lines, which are good tracers of stellar activity. 
While Ha is a pure absorption line in K2-3 and does not 
measurably vary, its Can lines, as for all M dwarfs, do have 
emission cor es. We quantify the C all H&K flux through 
the S index (Vaughan et al. 19781 and And that it varies 
by 30% on a timescale that is longer than than the ten-day 
period of Planet b. K2-3 has a lower average S index than 


j 1 0.6-day rotation period 
dBonfils et aL]|2013|), consistently with Prot > 10 days. We 
used our recent calibration (Astudillo-Defru et al. in prep.) 
to compute the index from S, and from the R'^x 

Prot in the same paper, we estimate Prot — 38 days. The 
K2 lightcurve varies with a < 2 mmag semi-amplitude and 
a ~40-day characteristic timescale, consistent with a fairly 
inactive star and with the estimated Prot- If that photo¬ 
metric variability is entirely due to a dark spot rotating in 
and out of view, the corresponding RV semi-amplitude is 
at most ~1.6 ms“^. 

In an effort to quantify activity-induced RV variations, 
we developed new activity metrics. We built an "active" 
and a "non-active" template from the third of the observed 
spectra with the highest and the lowest S-indices, respec¬ 
tively, computed two sets of active and non-active RVs, and 
used their difference ARV=RVact—RVnoact as an activity 
tracer. The ARV time series (Fig. varies by ±8 m/s for 
BJD—2 450 000 < 7070, and is consistent with zero after 
that date. We surmise that the early large variations result 
from a spot that not only occults part of the star but also 
imprints the overall spectrum with its spectral signature. 
That spot would have been present during the first ~ 25 
days of the observations and then disappeared. 


We computed our nominal RVs (Sect. 2.21 with a tem¬ 
plate that includes all spectra regardless of their S index, 
and it presumably has a sensitivity to activity that is inter¬ 
mediate between that of the active and non-active template. 
We account for its stellar activity sensitivity by introducing 
a term proportional to ARV in the radial velocity model, 
with the proportionality coefficient a adjusted as a free pa¬ 
rameter in the fit. 


3. Analysis 

We analyse this multi-planetary system with an n-body dy¬ 
namical model that describes the gravitational interactions 
between all components of the system and not just the 
pull of the central star. Unlike a Keplerian model, the dy¬ 
namical model therefore accounts for the photometric and 
radial-velocity changes induced by the mutual attraction of 
the planets, such as transit timing variations (TTVs), tran¬ 
sit duration variations (TDVs), or more generally, transit 
shape variations. Even though we detect no gravitational 
interactions, using a dynamical model helps constrain the 
masses and eccentricities by excluding values for those pa¬ 
rameters that would result in detectable interactions or in 
a highly unstable system. 

Our dynamic al model uses the well known MERCURY 
(Chambers 19991 n-body integrator code to compute the 
three-dimensional position and velocity of all system com¬ 
ponents as a function of time. The line-of-sight projection 
of the stellar velocity is compared to the HARPS radial ve¬ 
locities. The sky-projection of the planet-star separations 
are used, together with the planet-to-star-radius ratios and 
the limb darkening coeff icient, to model the K2 light curve 
(Mandel & Agol 20021. We couple the dynamical model 
with a Monte Carlo Markov Chain (MCMC) code, de¬ 
scribed in detail in Diaz et al. (2014||, to explore the poste¬ 
rior distribution of the paramieters. For each MCMC step, 
we run MERCURY three times with the same model param¬ 
eters: 
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- once over the 80 days of the K2 observations, using 
the Bulirsch-Stoer algorithm and a 0.02-day timestep 
to model the light curve with a 1 ppm maximum pho¬ 
tometric error 

- once over the (disjoint) 103 days of the HARPS obser¬ 
vations, again with the Bulirsch-Stoer algorithm and a 
0.02 days timestep, to model the radial velocities 

- once over a 1000-year interval, using the hybrid 
symplectic/Bulirsch-Stoer integrator with a 0.5-day 
timestep (l/20th of the shortest orbital period) 


The last run is not compared to any observation but used 
to ensure the short-term stability of the system, while re¬ 
jecting any MCMC step where two orbits intersect or where 
a planet comes within 0.05 au of the star. (Tidal forces are 
not included in MERCURY.) One would ideally like to evalu¬ 
ate stability on longer time scales, with 1000 years chosen as 
a compromise between ideality and computational expense. 

The physical parameters of the model are the stellar 
mass and radius, the coefficient of a linear limb-darkening 
law, the systemic velocity, the planetary masses, the planet- 
to-star radius ratios, and the plane tary orbital parameters 
(a, e, i, ui, n, and M; see Table A.21 at reference time 
2 456 812 BJD. To minimise correlations, however, we use 
a different parametrisation for MCMC jump parameters: 
orbital period (P), conjunction time of the first transit 
observed by K2 (Tg), radial velocity semi-amplitude (K), 
y/e cos uj, and y/e sin uj. Additionally, we fit a global light 
curve normalization factor, one multiplicative jitter param¬ 
eter for each data set and a, the am plitu de of the ARV 
term of the radial velocity model (Sec. 2.31. 

We use non-informative uniform priors for all MCMC 
model parameters except the stellar mass and radius, for 
which we ad opt the Gaussian distributions of [Crossfieid] 


et al. (20151 (we also adopt all parameters in common as 
the starting point of our chains) and the longitude of the 
ascending node for which we adopt a Gaussian distribu¬ 
tions with tr = 2° t o enforce the observed physical bias 
towards coplanarity (Fabrycky et al. 20141. We adopt a 
non-informative prior for the limb-darkening even though 
this widens the distributions of several parameters, because 
limb-darkening models have often been found to be inaccu¬ 
rate when the data quality is high enough to probe them. 
We ran 48 MCMC chains of 180 000 steps and combined 
their results as described in Diaz et al. (20141. 


4. Results 


Table |A)^ lists the mode and the 68.3% confidence inter¬ 
val for the system parameters. Figures [T] and |A^ show the 
radial velocity measurements and the transit light curves, 
together with their respective models. By fitting a stellar 
spectral energy distribution to the K2-3 photometry (Cross- 
field et al.|2015 Table 1) as described inlDiaz et al. (2014p, 
we obtain a 47.5 ± 6.0 pc distance and E(^b-v) < 0.056 
with 68.3% confidence. 


Our results agree with Crossfield et al. 


(20151 for all 


parameters in common. In addition to adding the radial 
velocity information, our analysis differs in that we neither 
impose a circular orbit nor use informative priors on i?*/a 
or the limb darkening coefficient. We also model the transits 
of all three planets jointly rather than one planet at a time, 
enforcing consistency in the derived stellar properties. Our 
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Fig. 1. Time series for, from top to bot tom: a) the HARPS 
radial velocities of K2-3, b) ARV (Sec. |2.3P , c) the RVs corrected 
from the aARV of the maximum likelihood model together with 
the dynamical model (the solid black line represents the median 
model, and the shades of grey represent the 68.3, 95.5, and 99.7% 
Bayesian confidence intervals), and d) the RV residual from the 
maximum likelihood dynamical model. 


approach measures the stellar density more accurately than 
can be inferred from the average spectrum. 

The times of the individual transit are measured with 
standard errors (which roughly translate to a Icr detection 
limits for TTVs) of 1.4, 2.8, and 2.9 minutes for planets b, c, 
and d, respectively. Over the short time span of the HARPS 
observations, the radial velocity of the dynamical model 
differs from that of a three-Keplerian model by at most 
27cms“^, which is well below the ~2ms“^ measurement 
noise. 

Owing to the low eccentricities and almost edge-on in¬ 
clinations, all three planets almost certainly undergo sec¬ 
ondary eclipses. Table [Ar2| lists the epochs and durations of 
these secondary eclipses. 

To evaluate the robustness of our mass measure¬ 
ments, we extended the radial velocity by adding a lin¬ 
ear drift, a fourth planet, a term proportional to the bi¬ 
sector of the spectra, and any combination of them. We 
also experimented with restricting the analysis to those 
radial velocities that seem unaffected by stellar activity 
(BJD—2 450 000 > 7070). They always produced similar 
velocity amplitudes (and masses) for Planet b, but a wide 
range of values for Planets c and d. We conclude that our 
data robustly measure the mass of Planet b, 8.4 ± 2.1 M 0 , 
but not those of c and d. Many more observations and very 
careful analysis will be needed to disentangle signatures of 
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Planets c and d from the activity signal. Fignrej^shows the 
RV signal from Planet b after removing the other terms of 
the model. 



Phase [P=10.05432 d] 


Fig. 2. HARPS radial velocities of K2-3 phased to the period 
of Planet b, with the Keplerian model (solid black line) overlaid, 
after removal of the Keplerians signals from Planets c and b and 
of aARV activity term of the maximum likelihood model. 


When inclnded in the model, the fonrth planet con¬ 
verges to a Pe ~ 100-day period and absorbs part of the 
residuals seen in Fig. H (bottom). Model comparison favours 
the four-planet model only marginally over the simpler 
three-planet model, with the Bayesian evidence estimators 


of Perrakis et al. (20141 andjChib & Jeliazkov (20011 giving 
odd ratios of 4.6±0.3 and 4.9±2.1, respectively, in favour 
of the four-planet model. More data will thus be needed to 
firmly establish whether additional planets orbit K2-3. 

Figure adds Planet b to the mass-radius diagram of 
the known small transiting planets. Planet b must contain 
rock with at most a 50% water envelope. 

Planets around M stars are ideal for characterisation 
follow-up, thanks to the favourable planet-to-star surface 
ratio. The other known low-mass planets transiting bright 


M stars have a gas envelope (GJ436b and GJ3470b, Butler 


et al. 2004 Bonfils et al. 


20121 or may even potentially 


consist of l66% water (GJ12140, Gharbonneau et al.|2009 1. 
Our measurement of the density of Planet b shows that it is 



either entirely rocky or mostly rocky with a water envelope. 
It is the first planet in this category, and the brightness of 
the star makes it a prime target for follow-up. 
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Appendix A: Additional figures and tables 


Table A.l. HARPS radial velocity measurements of K2-3. 



Fig. A.l. Periodogram of HARPS radial velocities of K2-3. 
The horizontal lines correspond to 1, 2, and 3a confidence inter¬ 
vals. The false alarm probabilities (FAP) of the main peaks are 
0.044% (20-30 days) and 0.0015% (8-10 days). The vertical lines 
indicate the period of Planets b, c, and d (from left to right). 
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57085.773295 

30.4726 

0.0025 

-2.0 

0.780 

57085.810298 

30.4797 

0.0024 

-2.6 

0.786 

57100.614169 

30.4874 

0.0033 

0.2 

0.774 

57100.822420 

30.4809 

0.0030 

-0.9 

0.822 

57101.660042 

30.4890 

0.0034 

-0.2 

0.576 

57104.607546 

30.4795 

0.0035 

-2.3 

0.932 

57117.543836 

30.4809 

0.0032 

2.3 

0.545 

57135.517340 

30.4782 

0.0024 

0.1 

0.838 

57135.673963 

30.4848 

0.0021 

1.2 

0.959 

57136.486899 

30.4828 

0.0024 

-1.7 

0.742 

57137.518119 

30.4849 

0.0032 

-0.3 

0.866 

57137.679348 

30.4819 

0.0027 

1.5 

0.659 

57138.523321 

30.4822 

0.0029 

-0.1 

0.867 

57138.676714 

30.4858 

0.0030 

0.4 

0.839 

57139.518338 

30.4765 

0.0029 

1.0 

0.835 

57139.679473 

30.4784 

0.0028 

0.3 

0.813 

57146.549907 

30.4749 

0.0020 

-0.8 

0.825 

57146.668220 

30.4745 

0.0034 

-1.4 

0.687 

57148.541552 

30.4691 

0.0030 

-0.1 

0.946 

57148.662503 

30.4717 

0.0033 

1.6 

0.619 


Notes. include a systematic error of 60 cms 
(Bonfils et al.|2013l. 
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Table A.2. Model parameters. Posterior mode and 68.3% credible intervals. The orbital elements have their origin at the star 
(Asteroidal parameters in the mercury code) and are osculating elements for reference time tref = 2 456 812 BJD. We have low 
confidence in the mass and density of planets c and d (Sec.[^. 


Parameter 

Mode and 68.3% credible interval 


Stellar mass, M* [Mq]* 

Stellar radius, [^gI* 

0.612 ± 0.086 

0.553 ± 0.041 



Stellar density, p* [p©] 

3.58 ± 0.61 



Surface gravity, logp [cgs] 

4.734 ± 0.062 



Linear limb darkening coefficient, u* 

0.573 ± 0.088 



Systemic velocity, 70 [kms“^]* 

30.48024 ± 0.00063 




Planet b 

Planet c 

Planet d 

Semi-major axis, a [AU] 

0.0775 ± 0.0039 

0.1405 ± 0.0067 

0.2086 ± 0.010 

Eccentricity, e 

[0, 0.12]* < 0.29* 

[ 0 , 0.08]* < 0 . 20 * 

(0, 0.09]* < 0.19* 

Inclination, i [°]*'^ 

89.59+°;24 

89.70 ± 0.20 

89.79 ± 0.15 

Argument of pericentre, lo [°] 

180 ± 70 

89 ± 110 

351 ± 66 

Longitude of the ascending node, n [°]* 

180 (fixed) 

180.0 ± 1.8 

180.0 ± 1.8 

Mean anomaly, M [°] 

230 ± 69 

21 ± 100 

338 ± 67 

Period, P [days]* 

10.05429 ± 0.00022 

24.6491 ± 0.0033 

44.5705 ± 0.0059 

Transit epoch. To [BJD-2,450,000]* 

6813.41817 ± 0.00082 

6812.2784 ± 0.0018 

6826.2272 ± 0.0027 

Transit duration, Ti 4 [hours] 

2.563 ± 0.039 

3.345 ± 0.079 

4.07 ± 0.12 

Scaled semi-major axis, a/Ri, 

30.74l°f 

55.901H" 

82.917 

Transit impact parameter, b 

0 

^■^^- 0.12 

0.3l7+t 

0.359l+°f 

Secondary eclipse impact parameter, bgE 

0.3l0l°;°'(® 

0 

'^•'^'^- 0.22 

'-'•'JU_o ,20 

Secondary eclipse epoch, [BJD-2,450,000] 

6818.407+0°“ 

6824.59 ± 0.85 

6848.72+^4 

Secondary eclipse duration [hours] 

2.55 ± 0.25 

q qy-|-0.40 
-0.17 

^•'^-‘-- 0.21 

Radial velocity semi-amplitude, K [ms“^]* 

3.60 ± 0.87 

^■'^^-0.37 

2.84 ± 0.90 

Radius ratio, Rp/R* 

0.03477 ± 0.00064 

0.02770 ± 0.00070 

0.02495 ± 0.00074 

y/eCOSLJ* 

n 9 qq"I" 0'33 
u.zc.o_o 065 

0.01 ± 0.19 

0.263lg;°f 

y/esinu* 

0 05"*”®'^^ 
^•^^-0.14 

0 . 1287 +° 

0.00 ± 0.16 

Planet mass, Mp [Mq ] 

8.4 ± 2.1 

2-i7:3 

11.1 ± 3.5 

Planet radius, i?p[R 0 ] 

9 078“^®'^^ 

Z.U/0_Q Qg3 

i.O4^_0 065 

1.53 ± 0.11 

Planet density, pp [g cm~^] 

‘*•'3^-0.76 

1 qo”I"3.0 

J--o^_0.96 

17.5 ± 6.3 

K2 long-cadence multiplicative jitter* 

1.020 ± 0.049 



HARPS multiplicative jitter* 

1.43 ± 0.13 



Amplitude of the ARV term, q* 

-0.21 ± 0.19 




Notes. 1*1 MCMC jump parameter, lb reflected with respect to i = 90°, the supplementary angle is equally probable. 1*1 68.3% 
highest density interval, lb upper limit, 99% confidence. Tq = t^ef ~ ^ (M — E + esmE) with E = 2arctan |^)^l^^tan [1 — iv)] |, 

P = K = 2 ^0 = 1-98842 xlO^o kg, Rq = 6.95508 xlO® m, M© = 5.9736 xlO^^ kg, R® = 6,378,137 m. 
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Fig. A.2. As in Figure [^ but for the K2 photometry of K2-3 (from top to bottom and left to right: planet b, planet c, and 
planet d). Each transit is centered relative to the linear ephemeris. 























































